Burgers Equation Revisited 



m 
o 
o 

(N 



V. Gurarie 

Department of Physics, Theoretical Physics, Oxford University, 1 Kehle Rd, Oxford 0X1 UP, United Kingdom 

(Dated: 14 March 2002) 

This paper studies the ID pressureless turbulence (the Burgers equation). It shows that rehable 
numerics in this problem is very easy to produce if one properly discretizes the Burgers equation. 
The numerics it presents confirms the 7/2 power law proposed for probability of observing large 
negative velocity gradients in this problem. It also suggests that the entire probability function for 
the velocity gradients could be universal, perhaps in some approximate sense. In particular, the 
probability that the velocity gradient is negative appears to be p « 0.21 ± 0.01 irrespective of the 
details of the random force. Finally, it speculates that the theory initially proposed by Polyakov, 
with a particular value of the "anomaly" parameter, may indeed be exact, at least as far as velocity 
gradients are concerned. 
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I. INTRODUCTION 

The problem of the randomly driven Burgers equation, 
or ID pressureless turbulence, attracted a lot of attention 
in the literature in the last decade 0,1111111111,1111. 
Indeed, it is a tantalizing problem. On the one hand, it 
incorporates many features one would expect from a real 
3D turbulent fluid. On the other hand, it is much sim- 
pler than the 3D turbulence of the incompressible fluid, 
usually studied in the framework of the Navier-Stokes 
equation. Its simplicity prompted some researchers to 
suggest that it will become the "Ising model" of turbu- 
lence. Yet many published papers later, we are still far 
from a complete solution to this problem. In fact, many 
of the aspects of the Burgers equation remain the sub- 
ject of an ongoing debate. Several competing methods 
exist in the literature, each with its own range of ap- 
plicability, which provide different, sometimes mutually 
exclusive answers to various aspects of the ID turbulence. 

The situation could be improved if a way existed to 
check various statements made about the Burgers equa- 
tion. Such a check could be provided by numerical sim- 
ulations. But the numerical simulations in turbulence 
are hard to conduct, they require large expenditure of 
resources, and the answers they provide are often too 
ambiguous to either conflrm or disprove a particular the- 
ory. The situation in Burgers turbulence seems to be 
no different. Despite several influential papers on the 
numerical aspects of Burgers turbulence, those studies 
failed to rule out completely or conflrm any of the com- 
peting theories currently on the market. The exception 
to this seems to be the paper of J. Bee where he 
confirmed the asymptotics of a certain probability dis- 
tribution function (discussed below) predicted earlier in 
Ref. 0. 

In this paper, I would like to point out that there ex- 
ists an alternative way to study the Burgers equation 
numerically. The method involves minimizing the Hopf- 
Cole functional instead of trying to solve a complicated 
nonlinear differential equation. That turned out to be 
extremely easy to do numerically. A ten line program. 



which can be written in half an hour, produces, after sev- 
eral minutes of computer time, clean graphs of probabil- 
ity distribution functions. With the help of this program, 
I report confirmation of some of the existing theories. I 
also hope that the method, being so simple, can be used 
by other researchers to do a quick check if their proposed 
theory indeed does not contradict the numerical experi- 
ment. 

The rest of this paper is organized as follows. In sec- 
tion II I formulate the problem of Burgers turbulence 
and summarize many of the known facts about it. In 
section III I show how Burgers equation must be dis- 
cretized in order to facilitate doing numerics. In sec- 
tion IV I discuss the statistics of the velocity gradients 
of the Burgers equation using the methods of Lagrangian 
trajectories. I also discuss the "anomaly" hypothesis of 
Polyakov and explain its meaning within the Lagrangian 
approach. In section V I discuss the numerical simula- 
tions, show that they are compatible with the asymptotic 
results of Refs. H H and show that they are also com- 
patible with the Polyakov hypothesis if the "anomaly" 
parameter (3 is chosen to be equal to 3/2. Finally in the 
last section VI I also show that the compatibility with 
the known results on the behavior of pinned charge den- 
sity waves 0,0,0,01 also forces the value /? = 3/2. I 
speculate that since the problem of pinned charge density 
waves and the asymptotics of Refs. relate to differ- 

ent properties of Burgers turbulence, there is a chance 
that the Polyakov's theory with (3 = 3/2 is indeed ex- 
act, at least as far as the calculation of probabilities of 
velocity gradients is concerned, unless a remarkable and 
unlikely coincidence is at play here. 



II. FORMULATION OF THE PROBLEM 

The Burgers turbulence problem can be formulated in 
just a few lines. Take a pressureless fluid in ID, introduce 
its velocity fleld u{x,t) where x is space and t is time and 
write down the Navier-Stokes equation 



dtU + u dxU — V d^u — f{x, t). 



(1) 
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Here v is the viscosity of the fluid which is taken to be 
very small, /(x, t) is a force acting on the fluid, which is 
taken to be random, white noise in time, and a smooth 
function in space 



(2) 



where F{x) is a function which is positive at x = and 
smoothly goes to zero when x ^ L, where L sets the 
force length scale. As far as boundary conditions are 
concerned, it is usually assumed that u is either periodic 
in space (in which case F has to be periodic as well), or 
the space is infinite and u at infinity. Note that the 
center of mass motion decouples in both cases, to give 



di 



dx u 



dx /, 



(3) 



which is but the standard Brownian motion. For this 
reason one usually selects the force / to be a deriva- 
tive of another function h{x,t), f = dxh. With this re- 
striction, the average velocity vanishes ( u ) = and it 
makes sense to define the root mean square velocity as 
Urms = y{u^- It measures the average amplitude of 
the fluctuation of the velocity field. A basic assumption 
in the theory of turbulence is that this quantity does not 
depend on the viscosity v as it is taken to zero. That al- 
lows to determine J7rms via simple hydrodynamic scaling 
to be f/rms oc (Li^(0))i 

Once the equation is written, one would like to cal- 
culate various averages and probabilities. A particular 
popular question focuses on the probability distribution 
function of the velocity at two different points. Define 
P(cr) as the probability of observing a velocity gradient 
a = dxu{x) in a particular point in space. Similarly 
r) is the probability of observing a difference in ve- 
locities V = u{x + r) — u{x) at two points separated by 
a distance r from each other. These quantities received 
a lot of attention in the literature, possibly because they 
are related to the turbulent advection of particles sus- 
pended in the fluid, which is a very important topic in 
real 3D turbulence. 

The following is a highly subjective historical overview 
of the progress made in understanding As always, 
the choice of highlights is very personal one, and I apol- 
ogize to the authors whose contributions are not men- 
tioned here. 

The first to address the behavior of Pd{v) were V. 
Yakhot and A. Chekhlov 0, who argued that it should 
not be symmetric under w — > —v. They were also the 
first to study this function numerically. 

Soon thereafter, A. Polyakov j2| suggested an approach 
to calculation of these functions. He suggested an equa- 
tion that Pd{v) should satisfy. The drawback of his ap- 
proach was that it was not just one equation but a family 
of equations parametrized by a parameter (referred to as 
"^-anomaly") whose value remained undetermined. At 
about the same time, Bouchaud, Mezard and Parisi Q 
attacked the problem using replica trick and the Gaus- 
sian variational ansatz. The problem with their approach 



was that it was not well suited to computing quantities 
such as P{a) and Pd{vY 

In subsequent work |j] of the author with A. Migdal, 
it was shown that at least the asymptotics of P{(7) and 
Pd{v) at cr, w — > -\-oo can be derived carefully in the saddle 
point approximation. It was shown that 



P(a) ocexp(-^), a»Gi«^ 

Pd(v) oc exp ( -g^^^ j , W > C^rms 



(4) 



where G = —dxF{x)\x=a- The drawback of the saddle 
point approach was that it allowed to determine the right 
hand side asymptotics only. But that asymptotics is ex- 
pected to depend on whether the force / of is Gaus- 
sian. The asymptotics is true only if the random 
force of the Burgers equation is Gaussian. But it does 
not have to be. For example, one could take the random 
force whose strength is always less than a certain value, 
/I < /max- Then the asymptotics Q would not be ex- 
pected to hold. But if it is so, then one could legitimately 
ask if at least some part of the functions P{v) and Pd are 
universal, that is, they are independent on the details of 
the random force /, except via the force length scale L. 
On top of it, since any random Gaussian force leads to 
Eq. 101), one could ask if at least for the Gaussian force, 
the entire functions P{<j) and Pd{v) are universal (that 
is, independent of the force correlation function _F(a;)). 

Some time later, an important paper by W. E, K. 
Khanin, A. Mazel and Y. Sinai 5] came out. It was 
argued in that paper that the asymptotics of P{a) at 
large negative a is 



Pia) oc 



1 



(7 2 



(5) 



The arguments leading to Eq. ||SJ) seemed to be well de- 
fined and valid independent of the details of the random 
force setup. Those arguments were later elaborated in 
a later paper by U. Frisch, J. Bee and B. Villone j^. 
In all of these papers, the asymptotics Eq. Q was not 
disputed. 

At the same time, these authors argued that the only 
universal feature of the probability distribution function 
P(o-) is its "left tail" Eq. The rest of the function was 
proposed to be nonuniversal, depending on the details of 
the random force. 

These claims were later checked numerically by R. 
Kraichnan and T. Gotoh (7(|, but they were not able to 
unambiguously confirm or reject the power law 10. Some 
time later J. Bee @ did a very thorough numerical anal- 
ysis of the "left tail" and concluded that the power law 
Eq. lO is indeed correct. 

Here I would like to report a numerical confirmation 
of the tail Eq. jSJ. However, I would also like to suggest 
that the left tail is not the only universal feature of the 
function P{<t). 
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III. DISCRETE BURGERS EQUATION 

To study the Burgers equation numerically, its space 
and time have to be discretized. One way to do that 
would be to write derivatives as differences. That is a 
very crude scheme, however. Burgers equation allows a 
much more gentle way of discretizing its time, which is 
based on the Hopf-Cole transformation. In this section I 
am going to present the scheme as it is, and show how it 
produces the Burgers equation in the continuum limit. 

Consider the energy function Ei(x) where x changes 
from to 27r and i is the integer index (which will later 
become time). Now write down the following definition 

Ei{x) = min [\ - cos {x ~ y) Ei^i{y)] + hi{x) . (6) 

0<y<27r 

Here hi is a random function, which I take as 

hi{x) = Ai cos {x - ipi) , (7) 

where Ai is a random variable (for example, random 
Gaussian with or taking values uniformly distributed 
from to some maximum Amax), and 0i is a random 
phase uniformly distributed from to 27r. The sign min 
in Eq. © should be understood as a minimization over 
y which varies from to '27r of the expression in the 
square brackets. Eq. © plays the role of a definition of 
Ei in terms of Ei^i. As a starting point, one can take 
Eq{x) = 0. Finally, the functions Ei{x) will automati- 
cally be periodic in x. 

The equation Eq. © is the discrete Burgers equation. 
To see that, let us take the continuum limit. In the nota- 
tions of ©, the continuum limit is achieved when hi <^ 1, 



Ei ^ 1. The minimization condition reads 

sm{x ^ y) - dyE,^i{y) ^ 0. (8) 

Since Ei is small, one can expand 

x^y + dyEi-i. (9) 

Substituting it back to lO and expanding Ei^i{x + y~x) 
in powers ol y — x, one finds 

E,{x) - E,^i{x) + i {d^E,^i{x) f = h,{x). (10) 

The final step is replacing the discrete index i with the 
continuous time t, to find 

dtE+^{d^Ef ^h, (11) 

which is the so-called KPZ equation. Introducing the 
velocity u ~ dxE, one gets the Burgers equation 

dtu + u dxU — dxh, (12) 

with dxh identified as the force /. In this equation, the 
viscosity term i^Uxx is absent. But it is easy to see that 
find a global minimum in © is equivalent to the infinites- 
imal viscosity term in (|12|) . following, for example, Feigel- 
man |^ . The derivation below is somewhat technical and 
is given for completeness only, so it is possible to skip the 
equations Eqs. H13|l - H15() without losing any essential in- 
formation. 

Introduce an alternative way of writing down © 



exp {—Ei{xi)) = Hrn J J^da^j exp 



j<i 



[1 — COS {xj — Xj^i) + hj{a 



(13) 



It is clear that the integral, in the limit of vanishing i' 
effectively minimizes the expression in the exponential, 
doing the work which is accomplished by the min sign 



in l^. Passing to the continuum limit, one finds the 
functional integral 



exp(-£;(y,r)) = lim^y Vx{t) expl-^J^ 



dt 



—x^ + h{x, t) 



, x{T) = y. 



(14) 



But the right hand side is none other than the propagator 
of the Schrodinger equation, which enables us to conclude 
that exp {—E) satisfies the Schrodinger equation and as 



a consequence, 

dtE + ^{dxEf 



vd^E — h, dtU + u dxU — vd^u — dxh, 

(15) 
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which concludes the derivation. (In the last line of (|15|l I 
replaced the arguments of the function E(y, T) by x and 
t.) 

All this suggests that the following steps must be taken 
to simulate the Burgers equation numerically. First, gen- 
erate the random "force" hi{x) using Ai can be 
taken, for example, uniformly distributed on the interval 
from to A,„ax- Then solve the relationship © recur- 
sively to find Ei{x). The actual value of Ei is not im- 
portant for the numerics of the Burgers equation. What 
is important is the value of y which minimizes the ex- 
pression in the square brackets of © for each value of x. 
Given the relationship between x and y, one finds for the 
Burgers velocity m = dxEi{x) 



Ui{x) = Ui_i{y) + d^hiix). 



(16) 



This is the discrete analog of the obvious equation 
-^u{x{t),t) = f{x,t) where ^x{t) = u{x{t),t), the so- 
called Lagrangian trajectory (see below). Finally, for the 
velocity gradients ai{x) — dlEi{x) one finds 



CTi(x) = cos {x - y) 



cos {x -y) + cri-i{y) 



-t- dlh,{x), 



(17) 

which is the discrete version of the Eq. H21(l considered 
below. The relations Eq. ©, Eq. and Eq. (QJJ are 
those one has to iterate to solve the Burgers equation nu- 
merically. To stay close to the continuum limit in time, 
one has to take Amax *C 1. And finally, to do the min- 
imization in © in practice one has to discretize space 
as well. With 2000 spacial discretization point, reliable 
data can be accumulated over about 10000 steps in time, 
which takes about 10 minutes computer time on a Sun 
workstation. The result of one such calculation, with 
^max = 0.01 is shown on Fig 1. 
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FIG. 1; A typical plot of u{x) at a fixed time t taken from an 
actual computer simulation 

The random force chosen according to fixes the 
function F to be F{x) oc cos(a:) and the characteristic 
length scale L is then of the order of the size of the sys- 
tem. If one wants to use a shorter range force, one could 
take 



hi{x) = A, 



oo 

E 

n— — oo 



exp ■ 



(x — 27rn -|- (j)iY 

l2 



(18) 



or any other similar function, as long as it is periodic in 
space X. 



IV. STATISTICS OF THE VELOCITY 
GRADIENTS 

The properties of the solution to Burgers equation has 
been discussed in many publication (see for example Q). 
Here I list those which are relevant for the present dis- 
cussion. 

In the Burgers equation, the viscosity v is supposed to 
be infinitesimally small. Yet it cannot be set to zero. The 
reason for that is, the solution to this equation develop 
shock waves, or sharp drops in the value of the velocity 
field. The width of this drops is of the order of l/v and 
they become discontinuities in u when v is taken to zero. 
In the middle of the drop fd^u has a finite limit as — > 0. 

To study the gradients of the velocity it is advanta- 
geous to introduce the notion of the lagrangian coordi- 
nate x(t). It is the point which moves together with the 
fluid 



dx{t) 
dt 



u{xit)). 



(19) 



A discrete version of the lagrangian trajectory was given 
by Eq. jnj. Once x{t) reaches a shock wave, it becomes 
"trapped" in it. For the purposes of this paper, it is 
natural to assume that x{t) ends as soon as it reaches 
the shock wave (this is precisely what Eq. does) . In- 
troduce the notion of the gradient of the velocity on a 
trajectory defined by x{t), 



(^{t) = dxu{x{t),t). 



(20) 



As a consequence of the Burgers equation, 

<7 + a^ ^dxfixit),t). (21) 

Notice that the v term can be neglected since, by con- 
struction, x{t) avoids shock waves and terminates once it 
reaches them. 

The equation Eq. (|21|1 can be interpreted as a Langevin 
equation with the random force g{t) — dxf{x{t),t). Be- 
fore proceeding further, one needs to establish the statis- 
tics of g{t). x{t) is correlated with / in a complicated 
way, and it is not a priori clear if those correlations can 
be neglected. 

Fortunately, the discrete Burgers equation discussed in 
the previous section provides an answer to this question. 
g{t) is none other than d^hi{x). But x is related to y via 
Eq. lO, and therefore, x knows nothing about the ran- 
dom phase (pi included in the definition of hi in Eq. |(7J). 
Therefore, g{t) is indeed just a white noise random force 



[g{t)git')) = GS{t^t') 



(22) 



where, as before, G = —d^F{x)\x=o ~ (Urms/L)^. Given 
a Langevin equation Eq. (|21|) . a Fokker-Planck equation 
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can be written down for the propability P(cr) 



d_ 



Gd_ 



dP 
'dt' 



(23) 



In the present context this equation was first written 
down in '10*1, but in that paper it was given an inter- 
pretation different from the one here. In fact, a shghtly 
more general form of this equation wiU be useful below. 
It has the form 



d_ 



G_d_ 



]P + (3P 



dP 
'dt 



(24) 



where j3 is an arbitrary parameter. At this stage, (3 = 
0. However, later a version of Eq. 1)24(1 with nonzero 13 
will be useful. Notice that the large a asymptotics of P 
obtained from the Eq. H24|) is Eq. Q , irrespective of the 
value of (3. 

Eq. (|21|l is unstable. As soon as a becomes large nega- 
tive, it quickly reaches infinity. This process is none other 
than the formation of a shock wave. A quick estimate of 
the probability P(cr) in this regime is easy to do. Once 
a is large enough and negative, one can neglect g{t) in 
Eq. 1)21(1 completely and solve the equation to find 



a{t) 



t-tn 



t < to 



The probability can be estimated to be 



P{a) (X dt S{a - a{t)) oc 



1 

7^' 



L 



(25) 



(26) 



This estimate was used in some publications to conclude 
that this should be the correct power law in the probabil- 
ity distribution function of the Burgers equation Eq. 
However, to say so would be premature. 

It is instructive to reproduce the power law Eq. 1(26(1 
with the help of the Fokker-Planck equation. To do that, 
observe that one has to solve this equation with an ab- 
sorbing boundary conditions at a — > — oo, corresponding 
to the trajectories disappearing when shock waves are 
formed. The standard substitution 



P 



exp 



(27) 



leads to the Schrodnger equation of a quantum mechani- 
cal particle moving in a potential U (a) = j (2G) — (/3 -|- 



2" 9^ 



a 

2G 



(/3 + l)aH/ = _ 



'dt' 



(28) 



The solution to this equation reads 

oo 

*(a) - ^ C„ exp i^Ent) ^'„(a), (29) 



where Cn are arbitrary coefficients. En are the eigen- 
states, and VP^ are the eigenfunctions of the Schrodinger 
equation satisfying 



G d^^n f cr 



2 9ct2 



2G 



~iP+l)a]^n^En^n- (30) 



At large times t ^ {Ei 
survives 



En 



only the ground state 



F(cr,t) oc e^p{-EQt) Po{a 
Po{<t) = exp 



3G 



t:>{Ei 
*o(ct)- 



(31) 



A detailed analysis of Eqs. 1(28(1 - 1(31(1 at arbitrary (3 can 
be found in the paper by S. Boldyrev flU. In the case of 
interest, f3 = 0. Ref. ^ij shows that in this case Eq > 0, 
and it confirms the asymptotics Eq. 1(26(1 for P(){a). The 
asymptotics at large positive a coincides with Eq. 10}. 
The fact that £^o > is not surprising. It means that the 
total probability / da P{cr) dicreases with time. This is 
due to elimination of trajectories which ended on shock 
waves. 

However, this cannot be the correct probability of ob- 
serving a velocity gradient cr at a given point in space. 
Indeed, introduce the density of trajectories p(x,t). It 
satisfies the continuity equation 



dtp + dx {up) = 0. 



(32) 



In particular, the density around the given Lagrangian 
trajectory x{t) is given by PL{t) = p(2;(i),i) and it satis- 
fies 



dt 



PL + o-pL = 



(33) 



It is clear that the probability that a particular trajec- 
tory happens to go through an observation point where 
a velocity gradient is measured is inversely proportional 
to the density of trajectories at this point. Thus a more 
physically relevant probability distribution function is de- 
fined as 



P[a) = { ^6{a-a{t)) 

PL 



(34) 



This function satisfies a Fokker-Planck equation Eq. 124(1 . 
but with (3 = \. In this case, it is possible to show that 
i?o < and 



P(cr) oc 



1 



r, cr-C 



C/r„ 



(35) 



n=0 



This asymptotics was proposed in the paper 0- But 
it cannot be the correct asymptotics either. Indeed, 
while the above analysis correctly removes the trajec- 
tories which end at the points where a new shock wave 
is formed, it fails to take into account those which fall 
into "mature" shock waves, the ones which formed in the 
past. That is why _Eo < and the probability grows with 
time. That comes from overcounting the trajectories. 
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Unfortunately, the Eq. H21() does not seem to contain 
enough information to determine if there is a shock wave 
nearby. If that is indeed the case, then the full Burgers 
equation Eq. ^ has to be solved. That is the core of 
the problem of Burgers turbulence, and the reason the 
solution to it seems so close and yet so hard to get. The 
full Burgers equation Eq. ^ leads to a full fledged 2D 
quantum field theory (via Martin-Siggia-Rose formalism, 
for example). That theory is presumably very hard, if not 
impossible, to solve. The "reduced" equation Eq. l(?T|l . on 
the other hand, leads to a quantum mechanics encoded in 
Eq. (|28|l . The solution to quantum mechanics problems 
are clearly within reach. The question is, therefore, if it is 
possible to extract all the relevant information about the 
Burgers equation from the Lagrangian equation Eq. (|21|l . 
At this stage of the present analysis, it does not seem 
possible. 

To move further, I will make the following leap of faith. 
The equation Eq. (|24ll was studied in Ref. llj for arbi- 
trary /3. It was found that while the large positive a 
asymptotics of P is always given by Eq. Q), the other 
asymptotics is given by 

P{a)^J^,a«-'^. (36) 

On the other hand, it is known from the analysis of 
Ref. Q and that the correct power of this asymptotics 
is 7/2. To match this asymptotics, I choose (3 = 3/2. 
Then I propose to solve the equation Eq. (|24l) and com- 
pare what it gives with the numerics. The results of the 
comparison is shown in next section. 

The equation Eq. (|24|l with an arbitrary (3 was pro- 
posed in 131 (but for Pd{v) as opposed to P{a)) where /3 
was referred to as 6-anomaly. But whether the present 
discussion has anything to do with the methods of Ref. 
is unclear. To derive Eq. H24|l with (3 — 3/2, one needs to 
show that the rate at which the trajectories fall into ma- 
ture shock waves is proportional to the extra ^crP term 
generated in the Fokker-Planck equation Eq. 124(1 . How 
to do that is not known to me. Next, the Eq. ((24|l can- 
not be exact. After all, it clearly breaks down even if 
the force is not Gaussian. However, it may be correct in 
some mean field sense, for moderate values of a. Finally, 
I'd like to note that for this value of (3, Eg < and the 
probability still grows with time. That must be because 
this description gives the relative probability of observing 
one value of a over another, not the absolute probabil- 
ity, and the overall normalization must be enforced by 
P P/U where U = J da P{a). 



V. NUMERICAL SIMULATIONS 

With the setup descibed in the preceeding sections, it 
is very easy to calculate P(cr) numerically. Unlike some 
of the previous work on this subject, I did not try to 
compute P{a) from the data. Instead, I computed the 



integrated probability density 

N{a) = f df, P(^). (37) 

Obviously its left tail is expected to go as N{a) oc 1/|(t| 4 . 
To plot it numerically, it is enough to collect the data into 
an array, to sort it in the increasing order, and plot the 
position of an array entry as a function of its value. Fig 
2. depicts N{(t) in the case when the force was chosen to 
be as in Eq. (jT)), Amax = 0.01, and J7rms « 0.02. 



N 




- 0.04 ^ -0.02 ^ ^ ^ ^ 0.02 ^ 0.04 ^ ' O 

FIG. 2: The integrated probability density N{a) as a function 
of a 

Fig 3. shows the left tail in the log-log format. For 
comparison, a straight line with the slope —5/2 is plotted. 
The agreement is striking. One could object that the 

LogA^ 



2 




FIG. 3: The plot of log A'^ as a function of log(— cr) for negative 
a. The straight line has a slope of —2.5 

power law in this graph does not extend over a sufhciently 
long range of a. To extend it further, one needs to do 
numerics at weaker force for longer periods of time (to 
keep close to the continuum limit even at large velocity 
fluctuations). It was not the task of the present work to 
do large scale numerical simulations, and I believe the 
agreement I have here is good enough (especially in view 
of the work reported in Ref. Q). 

Fig 4. shows the same probability function, but plot- 
ted together with the (numerical) solution to the equa- 
tion Eq. ^ with f3 chosen to be /3 = 3/2. To arrive at 
this picture, one parameter - the units of measure of cr - 
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N 



equation 




a 



FIG. 4: N{a) taken from the numerics (upper curve) and 
extracted from the solution to the Eq. 12411 with /3 = 3/2. 



had to be adjusted to make the graphs overlap as much 
as possible. The agreement is not perfect. Yet the curves 
are close enough to suggest that perhaps it is not coinci- 
dential either. One interesting number is the probability 
p that the velocity gradient is negative. It is especially 
interesting, because it is independent of the adjustments 
of the units of a. Numerically p « 0.21 ± 0.01. The 
solution to the Eq. (|24l) gives about 0.18. According to 
some claims in the literature, this number is not univer- 
sal. However, the value of p is reproducible whether I do 
the numerics with the long ranged force Eq. Q or shorter 
range force Eq. ()18|l with various values of L. So numeri- 
cally p appears to be universal, or perhaps approximately 
universal. 

To conclude, the numerical data does not confirm that 
the Eq. H24|) with the parameter /3 = 3/2 is the exact solu- 
tion to the problem. But it may be an approximate, per- 
haps in some mean field sense, solution to the problem, 
which reproduces the actual probability density rather 
close. I believe this question deserves further investiga- 
tions. 



VI. PINNED CHARGE DENSITY WAVES 

In the papers of the author with J.T. Chalker 0,^3 ^ 
remarkable correspondence was established between the 
randomly driven Burgers equation and the physics of 
pinned charge density waves. It would go beyond this pa- 
per to discuss charge density waves, but mathematically 
the correspondence can be formulated in the following 
way. Consider the functional 



S[x{t)] 



dt 



-iP' + h{x, t) 



(38) 



where h(x,t) is the same random function as the one in 
Eq. H13(l . Let us find a function xo{t) which, when substi- 
tuted for x{t) in Eq. (|38|l . gives an absolute minimum for 
the functional £. Obviously it satisfies the minimization 



^'xo dh{x,t) 



dx 



(39) 



x=xo{t) 



Now consider the energy cost of x{t) deviating from 
X{i{t). Writing x{t) k, XQ{t) -\- ip(t) we find the equation 
for the normal modes of oscillations about the absolute 
minimum, given by 



d^ ^ d^h{x,t) 
"dt^^ 



dx^ 



x=xo{t) 



Ip = U! tp- 



(40) 



With the problem thus set up, one needs to calculate the 
average number of modes with frequencies less than a 
given frequency w. Such a function is denoted as N{uj). 

This problem, which was first formulated in the con- 
text of charge density waves in Ref. received a 
substantial amount of attention in the literature (see 
Refs. H mill EE E IS- After initial disagreements, 
a consensus was built which gave for the function N{uj) 
the value 



N{lu) oc uj^ 



(41) 



for sufficiently small uj. This was first proposed in 
Ref. jl6l | and then the derivation was improved in sub- 
sequent publications until Ref. derived Eq. (|41|l in a 
systematic way. 

One of the results of Refs. 0,^3 relates the calcula- 
tion of N(uj) to the following specific question in Burgers 
turbulence. Consider a Burgers fluid Eq. |^ driven by a 
random force /(x,t) = dxh{x,t). Consider a Lagrangian 
trajectory a;o(i) which moves with a fluid and never gets 
absorbed by a shock wave for times t < T. Consider 
a velocity gradient a{t) defined in Eq. 120() . Let us find 
time intervals ti < t < t2 such that a{t) < within those 
intervals. Now let us calculate the probability that 



dt a{t) < log(w). 



(42) 



for some small value oicu. Then the function N{uj) of the 
charge density wave problem coincides with this proba- 
bility. 

Notice that knowing the tails of the probability distri- 
bution P{a) would not help calculating N{uj). This is 
because the typical functions a{t) which satisfy Eq. 
are not those for which a is very large negative, but rather 
moderate negative a{t) which however persist over lonj 
time intervals. Therefore, the results of W. E at al 
are useless if we were to calculate N{u!) with the help of 
Eq. ggi. 

Let us however use the methods of this paper, in par- 
ticular Eq. to calculate N{uj). To do so, we need 
to write down a Feynman path integral formalism which 
corresponds to Eq. H28|) and then use the Lagrange multi- 
plier method described in Ref. section VI C, subsec- 
tion 1 , corrected however for the presence of a parameter 
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/3 in Eq. HHJ- Then we find 

N{lu) oc a;2/3+2, (43) 

Now we know the exact answer to the problem, Eq. (|41|1 . 
If the methods described in this paper were to reproduce 
the exact answer, we have to choose f3 = 3/2. However, 
this is the same value of (3 required to reproduce the 
asymptotics of W. E et all 

Two possible explanations of this are possible. First is, 
perhaps the Eq. I|28|l with (3 = 3/2 indeed reproduces the 
right asymptotics of W. E et al and for some reason also 
reproduces the function N{uj) correctly, but this equation 
is still generally not correct and the probability distribu- 
tion derived with its help and shown on Fig. 4 is not 
correct either. Second is, the Eq. is indeed correct 
and gives the correct solution to the Burgers problem. 

The first explanation seems far fetched, as a remark- 
able coincidence must be at play to give rise to it. The 
second explanation sounds much more likely. And yet, 
no derivation of Eq. H28|) is known at this point. 

VII. CONCLUSION 

In this paper, I demonstrate that it is easy to obtain 
reliable numerical data on the behavior of the ID Burgers 



equation. I show that the data confirms the 7/2 tail of 
the probability distribution suggested in the literature. 
The data also suggests that the entire probability dis- 
tribution function could be universal, determined by the 
solution to the equation Eq. H24() with the "anomaly" 
parameter /3 = 3/2. I also discuss that it is possible to 
reproduce the known solution to the pinned charge den- 
sity wave problem by applying the techniques discussed 
here and choosing /3 = 3/2, which gives extra weight to 
these methods. 
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